Dipole-dipole interaction and polarization mode in BEC 
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We propose the construction of a set of quantum hydrodynamics equations for the Bose-Einstein con- 
densate (BEC) where atoms have electric dipole moment (EDM). The contribution of the dipole-dipole 
interactions (DDI) to the Euler equation is estimated. Quantum equations for the evolution of medium po- 
larization are constructed for the first time. The mathematical method we developed allows studying the 
effects of interactions on the evolution of polarization. The developed method may be applied to various 
physical systems in which dynamics is affected by DDI. A problem of elementary excitations in BEC, 
either affected or not affected by the uniform external electric field, is addressed using our method. We 
show that the evolution of polarization in BEC leads to the formation of a novel type of elementary exci- 
tations. Also, we consider the process of wave generation in polarized BEC by means of monoenergetic 
beam of neutral polarized particles. We compute the possibilities of the generation of Bogoliubov's modes 
and polarization modes by the dipole beam. 
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I. INTRODUCTION 

Since obtaining the Bose-Einstein condensate (BEC) in 
experiment with pairs of alkaline metal atoms, theoretical 
and experimental investigation of linear waves and nonlin- 
ear structures in it have been performed. The interest to the 
spinor BEC UJ-yD, the effect of magnetic moment on the 
BEC evolution yHg] and the influence of electrical polar- 
ization of atoms on the dynamic processes in BEC OLSHIl 
is rising in recent years. 

Many processes in quantum systems are determined 
by the dynamics and dispersion of elementary excitations 
(EE) 19|]. A law of dispersion of EE in degenerate dilute 
Bose gas was obtained by Bogoliubov in 1947 
Later, many authors studied the change of the Bogoliubov's 
mode which arise when interactions are more carefully in 
1 3 HI, a geometry of the system are complex Il4l[l5ll . In 
article 0] authors studied the influence of electric dipole 
moment (EDM) dynamics on dispersion of EE in BEC. The 
contribution of polarization in dispersion law of Bogoli- 
ubov's mode was obtained in [7]. The polarization waves 
in low dimensional and multy-layer systems of conduc- 
tors and dielectrics are considered in the papers [16l 171. 
Consequently, in BEC of atoms with EDM we expect the 
existence of polarizations wave; along with Bogoliubov's 
mode. It is particularly important to take into account 
the dynamics of EDM in investigations of exciton BEC in 
semiconductors. The method of QHD II 1811 was used in 
investigations of BEC lfl3ll and polarization waves in con- 
ductors and semiconductors lilal . Therefore, the method 
of QHD could be used to consider the possibility of arising 
polarization waves in BEC and calculating of the disper- 
sion law of these waves. 

For investigation of BEC, in system of particles with 
EDM or magnetic moments, various theoretical methods 
are used. The model equations, proposed for the descrip- 
tion of spinor and polarized BEC (such as the non-linear 
Schrodinger equation (NLSE) [Ql |3l ELlal ) are generaliza- 
tions of the Gross-Pitaevkii equation llllll . 



Other methods were applied for investigation of the po- 
larized BEC and other systems of particles, with EDM, 
along with NLSE. A hydrodynamic formulation of the 
Hartree-Fock theory for particles with significant EDM is 
considered in [19]. In paper [fl9ll the Euler-type equa- 
tion was derived, from the evolution of the density ma- 
trix. The EDM dynamics in dimer Mott insulators causes 
the rise of the low-frequency mode [20]. A two-body 
quantum problem for polarized molecules is analyzed in 
I 2lh . The existence of states with spontaneous interlay er 
coherence has been predicted in__n22fl in systems of po- 
lar molecules. Calculations, in I122ll . were based upon the 
secondary quantization approach, where the Hamiltonian 
accounts for the molecular rotation and dipole-dipole in- 
teractions. Superfluidity anisotropy of polarized fermion 
systems was shown in lT23ll and their thermodynamic and 
correlation properties were investigated n24n . The effect of 
EDM in a system of cooled neutral atoms which are used 
for quantum computing and quantum memory devices was 
analyzed in Il25ll . 

The characteristic property of BEC in a system of ex- 
citons inducing in semiconductors Il26ll it is a significant 
value of EDM of excitons. This leads to strong interac- 
tion of excitons with external electric field and emergence 
of the collective DDI in exciton systems. QHD method 
may be applied to such systems along with quantum kinet- 
ics based on the nonequilibrium Green functions [|27ll or 
density matrix equations Il28ll . 

Notable success has been reached in the Bose condensa- 
tion of dense gases |29|, 3(1. Consequently, we need the 
detailed account of short range interaction in investigations 
of dynamics of EDM. In this work we account the short 
range interaction (SRI) up to the third order on interaction 
radius (TOIR) [tl3ll . This approximation leads to nonlocal- 
ityofSRI EEE 

Electrically polarized BEC can interact with the beam 
of charged and polarized particles by means of charge - 
dipole and dipole-dipole interaction. Such interaction leads 
to transfer of energy from beam to medium and, conse- 
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quently, to generation of waves. In plasma physics the ef- 
fect of generation of waves by electron ||32|| or magnetized 
neutron ll33ll beam are well-known. In presented paper we 
consider similar effect in polarized BEC. 

We start from the equations of quantum hydrodynamics. 
System of QHD equations consists of continuity equation, 
momentum balance equation, eq. of polarization evolution 
and eq. of polarization current evolution. We made first 
principles derivation of this equation. For this purpose we 
used manyparticle Schrodinger equation. In this article we 
analytically calculate the dispersion properties of polarized 
BEC. We have shown that the dynamics of EDM leads to 
existence of new branch in dispersion law. Consequently, 
there are the waves of polarization in polarized BEC, along 
with Bogoliubov's mode. We obtain appropriate contribu- 
tion to the dispersion of Bogoliubov's mode related to DDI. 
Then, we consider the process of wave generation in po- 
larized BEC by means of monoenergetic beam of neutral 
polarized particles. We predict new method of generation 
of Bogoliubov's mode and a new type of modes, rest of the 
polarization waves. 

This paper is organized as follows. We introduce the 
model Hamiltonian in Sec. II, and present the momentum 
balance equation. Farther, in Sec. II, we derive the equa- 
tion of evolution of polarization and obtain the influence of 
the interactions on the evolution of polarization. In Sec. Ill 
we calculate the dispersion dependence of EE in BEC. The 
polarization evolution is taken into account. We obtain the 
contribution of polarization in dispersion of Bogoliubov's 
mode and show the existence of new wave solution. In 
Sec. IV we study the wave generation in polarized BEC 
by means of neutral polarized particle beam. In Sec. V we 
present the brief summary of our results. 



II. THE MODEL 

Explicit form of the hamiltonian of considering system 
in a quasi- static approximation is 



+ 



-y. 
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(i) 



The first term here is the operator for kinetic energy. The 
second term represents the interaction between the dipole 
moment d" and the external electrical field. The subse- 
quent terms represent short-range {7^ and dipole-dipole 
interactions between particles, respectively. The Green's 
function for dipole-dipole interaction is taken as G"f = 

V?Vf(l/r y ). 



For investigation of EE dynamic in polarized BEC we 
derive the system of QHD equations. This system of equa- 
tions consists of the continuity equation, the Euler's equa- 
tion and for the case of polarized particles the equations of 
polarization evolution and equation of field. The system of 
equations is derived by methods described in 111 311 . 

The first equation of a QHD equations system is the con- 
tinuity equation 



d t n(r, t) + d a (n(r, t)v a (r, t)) = 0. 



(2) 



The momentum balance equation for the polarized BEC 
has the form 

mn(r, t) (d t + vV)v a (r, t) + d p p a0 (r, t) 



4m 4m \ n(r, tj 



n{r,t)Td a n{r,t) + -T 2 cPAn 2 (r, t) 



pP(r,t)d a E p (r,t), 



where 



and 



T = y J dr(ry 



dUjr) 
dr 



T 2 = 



7T 



dr{r 



,dU(r) 



(3) 



(4) 



(5) 



30 J ' v ' dr 

In equation (O we defined a parameter T 2 as (01. This 
definition differs from the one in the work fl 1 3M - Terms 
proportional to h 2 appear as a result of usage of quantum 
kinematics. The first two members at the right side of the 
equation (0 are first terms of expansion of the quantum 
stress tensor. They occur because of taking into account of 
the SRI potential [7y. The interaction potential fTy deter- 
minates the macroscopic interaction constants T and T 2 . 
The last two members of the equation © describe force 
fields that affect the dipole moment in a unit of volume as 
the effect of the external electrical field and the field pro- 
duced by other dipoles, respectively. The last member is 
written using the self-consistent field approximation Ill8ll . 
p a/3 (r, t) is a tensor of the kinetic pressure, which depends 
on particles' thermal velocities and does not contribute into 
the BEC dynamics at temperatures near zero. 

Also, we have a field equation 



VE(r,t) = -47rVP(r,t). 



(6) 



In the case particles do not bear the dipole moment, the 
continuity equation and the momentum balance equation 
form a closed system of equations. When the dipole mo- 
ment is taken into account in a momentum balance equa- 
tion, a new physical value emerges, a polarization vector 
field P a (r, t). This causes system of equations to become 
incomplete. 
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Next equation we need for investigation of EE dispersion 
is the equation of polarization evolution 



d t P a {r,t) + d p R ap {r,t) = 0. 



(7) 



The equation (0 does not contain information about the 
effect of the interaction on the polarization evolution. The 
evolution equation of R al3 (r, t) can be constructed by anal- 
ogy with the above derived equations. Method of the equa- 
tions derivation is described in Appendix. Using a self- 
consistent field approximation of the dipole-dipole inter- 
action we obtain an equation for the polarization current 
R af) {r,t) evolution 



d t R af) (r,t) +d 7 ^R a/3 (r,i) 



-i? QT (r, t)v {r, t) - P a (r, t)v' s {r, t)w 7 (r, t) 



l h 2 

+ _d~<r al3 ~<(r,t) -d AP a (r,t) 

m 4m 2 



+ #_ &y ( d p P a (r,t)d in (r,t) | d,P a (r, t)d p n(r, t) 



8m 2 



n(r, t) 



n(r,t) 



-Td p (n(r,t)P a (r,t)) 



m 



m n(r,t) 



(8) 



Here r a/37 (r, t) represents the contribution of thermal 
movement of polarized particles into the dynamics of 
R a *(r,t). As we deal with BEC below, the contribution 
of r Q/37 (r, t) may be neglected. The last term of the for- 
mula ([8]) includes both external electrical field and a self- 
consistent field that particle dipoles create. Used in dSJ ap- 
proximations are described in Appendix. 

The first term in right side of Eq. (|8} describe the short 
range interaction. 

We can see various interactions are included in the equa- 
tions ([3| and (HJ additively. At a short distances among 
particles acts both SRI and dipole-dipole interaction. At 
large distances remain only dipole-dipole interaction. 

From Eq. ([8]l we can see that the change of polariza- 
tion arise from both the dipole-dipole interaction and the 
short range interaction. The SRI among particle leads to 
displacement of particles. Consequently, as particles has 
EDM, where are motion of EDM, e.i. changing of the 
R a ?(r,t). 

The terms that are proportional to H 2 are of quantum ori- 
gin as they are analogs of the Bohm quantum potential in 
the momentum balance equilibrium. 



III. ELEMENTARY EXCITATIONS IN THE POLARIZED 

BEC 

We can analyze the linear dynamics of elemental exci- 
tations in the polarized BEC using the QHD equations (0, 
©, ©, d7]l and dU). Let's assume the system is placed in an 
external electrical field E = E a e z . The values of concen- 
tration n and polarization P = kE for the system in an 
equilibrium state are constant and uniform and its velocity 
field v a (r, t) and tensor R al3 (r, t) values are zero. 

We consider the small perturbation of equilibrium state 
like 



n 



n + 5n, v a = + v c 



pa = pa + §po< jR af3 = Q + $ Ra ^ (Q) 

Substituting these relations into system of equations (O, 
©, (0, dSJ and (O and neglecting nonlinear terms, we ob- 
tain a system of linear homogeneous equations in partial 
derivatives with constant coefficients. Passing to the fol- 
lowing representation for small perturbations Sf 

5f = f(uj,k)exp(—zuj + ikr) 

yields the homogeneous system of algebraic equations. 
The electric field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the electric field, we come to the 
equation 

A • E x = 0, 



where 



2 h 2 k A Tk 2 n A P 2 k 2 
A = oo 2 - - — - H ; — - + Aira— — 



4m 2 



2m 



mn 



2vrTA: 4 P 2 



m 2 oo 2 — h 2 k 4 /4 + mTk 2 n — mT 2 k 4 n 
In this case, the dispersion equation is 
A = 0. 

Solving this equation with respect to oo 2 we obtain a fol- 
lowing results. 

The dispersion characteristic for EE in BEC can be ex- 
pressed in the form of 



" 2 = ^(-~rn ^ I 

2m \ 2 



h 2 k 4 

2m 



P 2 k 2 
n 



+ T 2 n r + 4vra- 




-Tn k 2 



T 2 n k i + 4ttct 



P 2 k 2 



8irTk A P 2 



(10) 



In contrast to the nonpolarized BEC, where only Bogoli 



nonpola: 

EE 



1 311 . a new wave solution ap- 



ubiv's mode exists 
pears in a polarized system due to the polarization dynam- 
ics. Bogoliubov's mode corresponds to a negative solu- 
tion of the equation (ITOl ). New wave solution it is a wave 
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of polarization. In general case presented by formula (fTOl 
the frequency of polarization wave in BEC depend on P 
and T. To investigate the BEC polarization effect on the 
Bogoliubov's mode and the dispersion characteristic of the 
new solution we analyze extreme cases of the formula (fTOb - 

Formula (ITOT ) demonstrates that taking account of the 
BEC polarization dynamics leads to a new solution. 

Let's start with the case when the effect of polarization is 
low compared to the contribution of short-range effects, i.e. 
of members proportional to T and Y 2 . If so, then formula 
dTUl) takes the form 



2 1,4 



LJ, 



h 2 k 4 Tn k 2 T 2 n k 



4m? 



+ 



771 



77? 



4irP 2 k 2 T(<7 



1) - 2oT 2 k 2 



mn 



to. 



h 2 k 4 

4th 2 



+ 



4nP 2 k 2 T(a 



T - 2T 2 k 2 

1 Tn k 2 

2 771 

1) - 2aT 2 k 2 



(11) 



mn 



T - 2T 9 k 2 



(12) 



Here we use indexes "B" for Bogoliubov's mode and "P" 
for new polarization mode. When deriving formulae (fTTt 
and (fT2l) we expanded a sub-radical expression in ([Tot and 
took only first two terms of the expansion to estimate the 
influence of polarization on the wave dispersion. 

As it follows from equations © and ©, changes in po- 
larization can occur due dipole-dipole interactions as well 
as to SRI and to quantum Bohm potential, i.e. to members 
proportional to h 2 . That's the reason for existence of other 
ways of polarization change when the contribution of equi- 
librium polarization P is neglectable. In the latter case, 
taken in a linear approximation, changes in the polariza- 
tion do not affect the concentration evolution. 

Formulae ([Tol are valid even in the absence of the ex- 
ternal electrical field when equilibrium polarization equals 
zero P = 0. Equations ([Toll in that case take the form 



h 2 k 4 

4m 



Wb = — I 



tO, 



Tn k 2 + T 2 n k 4 




^~ Tti £; 2 



(13) 



(14) 



The waves of polarization could be existed at the absence 
of external electric field E = 0. In this case the equilib- 
rium state is no polarized and dipole direction of particles 
is distributed accidentally. 

If the contribution of equilibrium polarization into BEC 
dispersion ([Tol is comparable to the contribution of SRI 
in the third order of the interaction radius T 2 relationships 
dT0l l are transformed into 



2 h 2 k 4 Tn k 2 T 2 n k 

4m z 771 771 



2 1,2 



+ 4tt : 



mn 



(15) 



K z k 
4m 2 



lTn k 2 A P 2 k 2 . 
h 4ir (cr 

2 771 77171 



(16) 



Using Feshbach's resonance 03411 we can transform the 
SRI potential in such a way that T = while T 2 7^ 0. 
Formulas ( fTOl in this situation turn into 



U), 



h 2 k 4 T 2 n k 4 
+ ■ 



4777 1 



in 



h 2 k 4 A P 2 k 2 
+ 4ira— — . 



4m 2 



mn Q 



(17) 



(18) 



In this paper we primarily focus on the influence of BEC 
polarization on its dispersion characteristics. So, let's con- 
sider the case when the contribution into the dispersion of 
the SRI at the first order of the interaction radius, i.e. terms 
proportional to T, is comparable to the contribution of po- 
larization, and their total effect is much greater than the 
contribution of terms proportional to T 2 . Here, ( fTOl turns 
into 



h 2 k 4 Tn k 2 ({a- 1) 



4m 2 



777 



(7 



(19) 



to. 



h 2 k 4 A P 2 k 2 

+ 4tt<j— — 



4t?i 2 



717710 



l T7iofc 2 / (a + 4) 
2 777 V 2a 



(20) 



The the first order of the interaction radius interaction 
constant for dilute gases has the form 



T 



4nh 2 a 



777 



where a is the scattering length SL llllLll3ll . The value T 2 
may be expressed approximately as 

0a 2 T 



where is a constant positive value about 1 , which depends 
on the interatomic interaction potential. Finally, T 2 takes 
the form 



To 



7t8h 2 a 3 
2m 



(21) 
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IV. GENERATION OF WAVES IN POLARIZED BEC 

In that section we consider the process of wave genera- 
tion in BEC by means of the beam of neutral polarized par- 
ticles. The interaction between beam and BEC has dipole- 
dipole origin. 

To get the dispersion solution we use the system of QHD 
equations for each sort of particles (0, (0, (0, © and the 
equation of field ©. The equilibrium state of system is 
characterized by following values of the BEC parameters: 



n 



n Q + Sn, v a = + v a , 



P a = p** + 5P a ,R ap = + 5R ap 
and values of the beam parameters: 

rib = n ob + 5n b , v b = U5 za + 5v b , 



(22) 



pa — pa < X pa 



R, 



a/3 



+ SB? 



(23) 



The polarization P£ is proportional to external elec- 
tric field E%. We consider the case then E n = 



[E sinip, 0, E cos(p\. In this case the tensor R^ b has 
only two unequal to zero elements: R^ b = R ob simp and 
RH = R ob cosip. For the process, under consideration the 
dispersion relation is: 



1 + 4vr/c 2 



p2 



id* 



h 2 ^ i r 

4m 2 1 2m 



fc 2 



h 2 k 4 Tn k 2 T 2 n k i 
4m 2 



and 



in 



, 2 _^k 2 

Di 



111 



(27) 



(28) 



where i is the index of sorts of particles, the BEC or the 
beam. 

The equation (l25l l has two beam related solutions, in the 
absence of BEC medium: 



h 2 k 4 

u = k z U ± J + a b u 2 Db . 



(29) 



We will consider the possibilities of instabilities for the 
case of low-density beam, the limit case uj Db ~ n ob — > 0. 
In this case we can neglect the last term in square root in 
( [291 . The resonance interaction beam with the BEC realiz- 
ing if 

hk 2 

k z U±—=uj(k), (30) 
2m b 

and could lead to instabilities. The quantity uj(k) is the 
dispersion of BEC modes (TTOl . The frequency in this case 
can be presented in the form 



I TT J_ ^ S 

uj = k z U ± h dui. 

2m b 



Lets consider two limit cases. 



(31) 



Tk 2 /(2m 2 



a \ 



h 2 k 4 , Tn k 2 
4m 2 m 



T 2 n fc 4 



mn 0/ 



small frequency shift limit 



In the limit case 



+ 



(w - k z uy 



fe 2 fc 4 

4m 2 



2(w - k z U)Pobh(P ob U - R ob ) 



m b n b[ (w - k z U) 2 



Using relation P bU — R ob 
tion (l24l and obtain 



4m? 



=0. 



(24) 

0, we can simplify the equa- 



1 + 



u 2 D ( Tn k 2 /(2m) 



a 







k*U) 2 - 

In this formula the following designations are used 

Tn Q k 2 



4m 2 



2m 



(25) 



(26) 



5uj -C hk 2 /m 
the frequency shift obtain in the form: 



5u 2 



2o b m b m 2 uj 2 Db (uj 2 — uj 2 ) 2 (uj 2 — u 2 ) 2 
' ojoj 2 D hT 2 n 2 k 6 W 



where 



W = 2u 2 



2ma 



2^2 



(32) 



(33) 



(34) 



Tn k 2 

and the frequency u determined with formula (TTOt . The 
instabilities take place in the case buS 2 < 0. The sign of 
8u 2 is depend on the sign of W. 

For the case resonance interaction of beam with the 
waves in BEC there are instabilities, for the first beam 
mode in rt29l at W < and for the second beam mode 
in formula (l29l if W > 0. For the polarization mode W is 
positive. It means that the interaction of polarization mode 
with the second beam related mode results in the instabil- 
ity. For the Bogoliubov's mode the sign of W depend on 
a. 

We can consider the following cases: 
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(i) the contribution of equilibrium polarization to u(k) 
is dominant; then W > 0; 

(ii) the dominant contribution to cj(fc) results from the 
term in u){k) which is proportional to T, and, equilibrium 
polarization and SPI in TOIR give comparable contribution 
tou)(k). In this case there is 



0"o 



1 + 



T 2 n fe 4 /m 



Tn k 2 /(2m) + 2T 2 n k i /m - 4uj 



. The sign of W varies at a = a a . The dependence of W 
from a is presented in the table 1 . 



large frequency shift limit 



In the limit case 



we have 



5uj » hk 2 /m 



(35) 



e, , t 3/ (T b m 2 UJ 2 Db (uj 2 - Ujf) 2 (uJ 2 - u\) 2 

dUJ = £\ 2/41 ixr I ' ( 36 ) 

where £ equal to £i = \/l for W > or £_i = for 
W < 0. Evident form of quantities ^ and £_i is 



V. CONCLUSION 

In this work we developed the method for description 
of dynamics of polarized BEC. This method accounts for 
the effect of polarization on changes in the concentration 
and in velocity field, which are determined in general by 
the continuity equation and Euler's equation. We derived 
the evolution equations of the polarization and the polar- 
ization current. The equations derived contain information 
about the influence of the interactions on the polarization 
evolution. We studied the effect of polarization on the BEC 
dynamics and the influence of SRI on the polarization evo- 
lution. An expression of SRI contribution in the equation 
of polarization current evolution via concentration, polar- 
ization and the SRI potential T = T(£/y) was derived. 
With the assumption that the state of polarized Bose parti- 
cles in the form of BEC can be described with some single- 
particle wave function. Changes in polarization due to SRI 
are shown to be determined at the first order of the interac- 
tion radius by the same interaction constant that occurs in 
Euler's equation and Gross-Pitaevskii equation. 

The dispersion of EE in the polarized BEC was analyzed. 
In the article, we show that polarization evolution in BEC 
causes a novel type of waves in BEC. The effect of polar- 
ization on the dispersion of the Bogoliubov's mode and the 
dispersion of a new wave mode were studied. 

We show the possibility of wave generation in polarized 
BEC by means of monoenergetic beam of neutral polarized 
particles. 



1, 



and 



6 = [1, 



1 + 1-*V3, 



The considerations concerning the sign of W, reflected 
in the table 1, are valid also for the limit condition ( |35l l. 

From the formulas ( 133) and ( |36l l we can see that external 
neutral particles beam leads to instabilities for both Bogoli- 
ubovs and polarization waves. 



TABLE I. In this table the sign of the W are presented for the 
Bogoliubov's mode when SRI is prevailed. 





T > 


T < 




(attraction) 


(repulsion) 


<J> <7 


+ 




a < a 




+ 



APPENDIX 

At derivation of system of QHD equations we need to 
write evident form of Hamiltonian of dipole-dipole in- 
teraction. In some works the dynamics of the magnetic 
dipole moment and of the EDM J3i l35ll are analyzed 
in similar ways. Usual expressions of a Hamiltonian for 
dipole-dipole interaction are equal for electric and mag- 
netic dipoles: 



H, 



dd 



However it has been shown by Breit [36] that a Hamil- 
tonian for spin-spin interaction, and, as a consequence, for 
the interaction of magnetic moments contains a term that 
is proportional to Dirac (5-function 8{r 1 — r 9 )d"d % T he 



coefficient of the <5-function has been refined later Q7II so 
that the Hamiltonian is in accord with Maxwell's free equa- 
tions, such as divB = 0. The resultant expression for the 
spin-spin interaction Hamiltonian is: 



H. 



IXfX 



47r5 a/3 (5(r 12 )+V?V?(l/r 12 ) U?/^. 
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Thus, the differences in the dipole-dipole interactions of 
electric dipoles, and magnetic dipoles, must be taken into 
account in the development of theoretical field apparatus. 

The Schrodinger equation defines wave function in a 
3N-dimensional configuration space. Physical processes 
in systems that involve large number of bodies occur in 
a three-dimensional physical space 113 811 . That's why a 
problem evolves of obtaining a quantum-mechanical de- 
scription of a system of particles in terms of material fields 
e.g. concentration, momentum density, energy density and 
other fields of various tensor dimension that are defined in 
a three-dimension space. 

In equations (O, © polarization occurs in the form of 



P a (r,t) 



J dR^^ 



- ri W(R,t)d?il>(R,t), (37) 



where r< is the coordinate operator for i-th particle, dR 



Up=i dr p . 

In right side of equation ([8| where is the force-like field 
F al3 (r, t) which give rise to evolution of the polarization 
current R a "(r, t). In general form, for F a "(r, t), we can 
write: 

F ap (r,t) = -—d^ a0 Vr,t) +—D a Vr,t)d p EVr,t). 
m m 

(38) 

A tensor 

D*?(r,t) = J dR S(r - r,)(^df t)if>(R, t), 

(39) 

occurs in the term that represents a dipole-dipole interac- 
tion and an interaction of the dipole with external electrical 
field. Based on the reasons of dimensions this value can be 
approximately presented as 



D a0 (r,t) = a 



n(r,t) 



(40) 



A SRI causes the tensor S"^ 7 (r, t) to occur in the equa- 
tion dSJ. Taken at the first order of the interaction radius it 
has the form 

Z^(v,t)=~JdRj26(r-*ii) 



4 r l Wity t)d?t/>(R, t) . (41) 



dr. 



Tensor X a/37 (r, t) describe the influence of SRI on evolu- 
tion of polarization. Eq. (|4H describe the SRI. 

If we apply the procedure described in ffl~3Tl to the cal- 
culation of the quantum stress tensor, and neglect the con- 
tribution of thermal excitations, £ Q/37 (r, t) for BEC takes 
a form of 

1 



Formula (l42l is obtained for the case where particles lo- 
cated in state with the lowest energy, which could be de- 
scribed by one particle wave function. This state may be 
the product of strong interaction. 

Tensor £ Q/37 (r,i) is, therefore, like the quantum stress 
tensor <r Q ^(r, t) in the momentum balance equation (O, 
dependent on T at the first order of the interaction radius 
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^BEC 



(r,i) 



T6^n(r,t)P a (r,t). (42) 
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